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Abstract — We extend the notion of memristive systems to ca- 
pacitive and inductive elements, namely capacitors and inductors 
whose properties depend on the state and history of the system. 
All these elements show pinched hysteretic loops in the two 
constitutive variables that define them: current-voltage for the 
memristor, charge-voltage for the memcapacitor, and current-flux 
for the meminductor. We argue that these devices are common at 
the nanoscale where the dynamical properties of electrons and 
ions are likely to depend on the history of the system, at least 
within certain time scales. These elements and their combination 
in circuits open up new functionalities in electronics and they are 
likely to find applications in neuromorphic devices to simulate 
learning, adaptive and spontaneous behavior. 

Index Terms — Memory, Resistance, Capacitance, Inductance, 
Dynamic response, Hysteresis. 



I. Introduction 

CIRCUITS elements that store information without the 
need of a power source would represent a paradigm 
change in electronics, allowing for low-power computation 
and storage. In addition, if that information spans a continuous 
range of values analog computation may replace the present 
digital one. Such a concept is also likely to be at the origin of 
the workings of the human brain and possibly of many other 
mechanisms in living organisms so that such circuit elements 
may help us understand adaptive and spontaneous behavior, 
or even learning. 

One such circuit element is the memory-resistor (memristor 
for short) which was postulated by Chua in 1971 by analyzing 
mathematical relations between pairs of fundamental circuit 
variables [1|. The memristor is characterized by a relation 
between the charge and the flux, defined mathematically as the 
time integral of the voltage, which need not have a magnetic 
flux interpretation. This relation can be generalized to include 
any class of two-terminal devices (which are called memristive 
systems) whose resistance depends on the internal state of the 
system [2|. 

Many systems belong to this class, including the thermis- 
tor (3) (whose internal state depends on the temperature), 
molecules whose resistance changes according to their atomic 
configuration [4|, or spintronic devices whose resistance varies 
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Fig. 1. (Color online) Symbols of the three different devices defined in this 
work: memristor, memcapacitor, and meminductor. In general, these devices 
are asymmetric and we use the following convention when possible: when a 
positive voltage is applied to the second (upper) terminal with respect to the 
terminal denoted by the black thick line, the memory device goes into a state 
of high resistance, capacitance or inductance, respectively. 



according to their spin polarization (5), 0. Recently, mem- 
ristive behavior and memory storage have been reported in 
solid-state TiC<2 thin films (7), O, where the change in 
resistance is realized by the ionic motion of oxygen vacancies 
activated by current flow. Additionally, memristive behavior 
has been demonstrated in VO2 thin films, where the memory 
mechanism is related to the insulator-to-metal transition in 
these structures J9). Finally, memristive behavior has been 
identified by two of the present authors (MDV and YVP) as 
a possible mechanism in the adaptive behavior of unicellular 
organisms such as amoebas fit)! . 

All these examples show the ubiquitous nature of mem- 
ristive systems. In fact, it should not come as a surprise that 
many of the above examples refer to nanoscale systems, whose 
resistance is likely to depend on their state and dynamical 
history, at least within (possibly very short) times scales 
dictated by the fundamental state variables that control their 
operation Q. 

In this paper we show that the above concept of memory 
device is not necessarily limited to resistances but can in fact 
be generalized to capacitative and inductive systems. Quite 
generally, if x denotes a set of n state variables describing the 
internal state of the system, u(t) and y(t) are any two com- 
plementary constitutive variables IfTTIl (i.e., current, charge, 
voltage, or flux) denoting input and output of the system, and 
g is a generalized response, we can define a general class of 
nth-order u-controlled memory devices as those described by 
the following relations 



y(t) 



g(x,u,t)u(t) 
f(x,u,t) 



(1) 
(2) 
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where / is a continuous n-dimensional vector function, and 
we assume on physical grounds that, given an initial state 
u(t = to) at time to, Eq. admits a unique solution |[T2l . 

Memcapacitive and meminductive systems are special cases 
of Eqs. (Q~|i and 0, where the two constitutive variables that 
define them are charge and voltage for the memcapacitance, 
and current and flux for the meminductance (see Fig. Q] for 
the symbols we use for these devices). In this paper, we 
will discuss several properties of these systems which make 
them unique compared to their "standard" definitions and 
will provide few examples. Combined with the general class 
of memristive systems and/or standard circuit elements these 
devices may lead to new functionalities in electronics. 

II. Memristive systems 

For completeness, let us introduce first the notion of 
memory -resistive systems 0. We also provide a simple 
analytical example of memristive behavior which has been 
employed in describing "learning circuits" iflOl . For a broad 
range of memristive system models, see, e.g., Ref. (14*1 . 

Definition - From Eqs. and 0, an nth-order current- 
controlled memristive system is described by the equations 

V M (t) = R(x,I,t)I(t) (3) 
x = f(x,I,t) (4) 

with x a vector representing n internal state variables, Vm (t) 
and lit) denote the voltage and current across the device, and 
R is a scalar, called the memristance (for memory resistance) 
with the physical units of Ohm. The equation for a charge- 
controlled memristor is a particular case of Eqs. (0 and (0, 
when R depends only on charge, namely 

V M =R(q(t))I, (5) 

with the charge related to the current via time derivative: / = 
dq/dt. Note that the analytical equation describing the Ti02 
device derived in the work by the Hewlett-Packard group [7 1 
has precisely this form, and it therefore represents an ideal 
memristor. 

We can also define an nth-order voZtage-controlled memris- 
tive system from the relations 

J(t) = G{x,V M ,t)V M {t) (6) 
x = f(x,V M ,t) (7) 

where we call G the memductance (for memory conductance). 

Properties - Several properties of memristive systems were 
identified by Chua and Kang in their seminal paper [2) . Here, 
we highlight just a few of them. Considering devices with 
R (x, I, t) > in Eq. 0, it has been proven that these devices 
are passive 0. Another important feature, following from 
the above inequality, is the so-called "no energy discharge 
property" which is related to the fact that a memristive system 
can not store energy, like a capacitor or an inductor. As a 
manifestation of these characteristics, we notice from Eq. 
[or Eq. (0] that Vm = whenever I — (and vice versa). In 



addition, for a periodic current input, the memristive system 
shows a "pinched hysteretic loop" 0. 

Moreover, a memristive system behaves as a linear resistor 
in the limit of infinite frequency and as a non-linear resistor 
in the limit of zero frequency provided that x = f (x , I) = 
in Eq. (0 (respectively x — f (x,Vm) = in Eq. 0) 
has a steady-state solution. These two properties are easy to 
understand on physical grounds. Irrespective of the underlying 
physical mechanisms that define the state of the system, at very 
low frequencies, the system has enough time to adjust its value 
of resistance to a momentary value of the control parameter 
(either current or voltage), so that the device behaves as a 
non-linear resistor (from Eq. or Eq. (0). On the other 
hand, at very high frequency, there is not enough time for any 
kind of resistance change during a period of oscillations of 
the control parameter, so that the device operates as a usual 
(linear) resistor. 

Finally, it is worth pointing out that, when a current flows 
in a resistor, at the microscopic level local resistivity dipoles 
form due to scattering of carriers at interfaces lfl5l . Namely, 
charges accumulate on one side of the resistor with consequent 
depletion on the other side, with possible accumulation and 
depletion over the whole spatial extent of the resistor fl6l . 
These dipoles develop dynamically over length scales of the 
order of the screening length [17], namely within length 
scales over which electrons can effectively screen the electric 
field of a local charge. In other words, the self-consistent 
formation of local resistivity dipoles takes some time and is 
generally accompanied by some energy storage related to the 
capacitance of the system. We thus expect that in nanoscale 
systems memristive behavior is always accompanied to some 
extent by a (possibly very small) capacitative behavior and 
as such the hysteresis loop may not exactly cross the origin 
at frequencies comparable to the inverse characteristic time 
of charge equilibration processes. This is, for instance, found 
in the transverse voltage response of the spin Hall effect in 
inhomogeneous semiconductor systems |[T8l . 

Examples - Several systems have been found to satisfy the 
above properties 0, 0, 0, 0, 0, 0, 0, and many 
more are likely to be found which fit this classification. In 
particular, memristive behavior is a property of thermistors 
0, molecular systems [4], spintronic devices [5] and thin 
film nanostructures [7|, [8|, [9|, as well as numerous examples 
presented in Ref. Ifl4l . 

As an example, let us consider a model of a voltage- 
controlled memristive system used by two of the present 
authors (MDV and YVP) in Ref. (TT5). In this model, the 
rate of system resistance change (R changes between two 
limiting values i?i and R2) is characterized by a constant a 
when \Vm\ < Vr and (3 when \Vm\ > Vr, where Vm is the 
voltage applied to the memristor and Vr is a threshold voltage 
[19|. Mathematically, the response of this memristive system 
is described by Eqs. (0 with 

G = x~ l , (8) 
x = -(PVM + 0.5(a-/3)[\VM + VT\-\VM-Vr\}) 
x [9 (V M ) 9{x-R{) + (-V M ) 6 (R 2 - x)] , (9) 
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Fig. 2. (Color online) Simulation of a voltage-controlled memristive system. 
This plot was obtained by a numerical solution of Eqs. \619\ for the case of a 
memristive system M connected directly to an alternating (ac) voltage source 
V (t) = Vn sin (2nujt) as shown in the inset using the following parameter 
values: K a = 20Q, R 2 = 500H, V T = 0.5V, a = -500V/(Q-s), = 2a, 
Vq = IV. In our model, the rate of resistance change is defined by the 
function /(Vm) = - (0V M + 0.5 (a - 0) [\V M + Vr\ - | V M - Vr\]) 
entering Eq. (9) and schematically shown in the inset for positive values 
of a and 0. At \Vm\ < Vt me slope of /(Vm) is determined by a and at 
\Vm \ > Vt it is determined by 0. 




Fig. 3. (Color online) Schematics of a pinched hysteresis loop of a memca- 
pacitive system. The energy added to/removed from the system (f Vc(q)dq) 
is the area between the curve and the q axis. The areas of shaded regions 
U\ and U2 give the amount of added/removed energy in each half-period. 
The signs of U\ and U2 are determined by the direction on the loop. For the 
direction shown here, U2 is positive and U± is negative. 



and a charge-controlled memcapacitor when Eqs. ( fl2t and ( TPS! 
can be written as 



where #(■) is the step function. The simulation results 
depicted in Fig. [2] clearly show a pinched hysteresis loop 
and a hysteresis collapse with increasing frequency uj of the 
alternating (ac) voltage source V(t) = Vq sin (2iru)t). 

We can now extend the above definitions to capacitances 
and inductances. The memory devices that result share many 
characteristics of memristive systems, but with a fundamental 
difference: they store energy. 

III. Memcapacitative systems 

Definition - We define an nth-order voltage-controlled mem- 
capacitive system by the equations 

q(t) = C(x,V c ,t)V c (t) (10) 
x = f(x,V C) t) (11) 

where q(i) is the charge on the capacitor at time t, Vcif) 
is the corresponding voltage, and C is the memcapacitance 
(for memory capacitance) which depends on the state of 
the system. Similarly, we can define an nth-order charge- 
controlled memcapacitive system from the equations 

V c (t) = C^&qrfqit) (12) 
x = f(x,q,t) (13) 

where C _1 is an inverse memcapacitance. 

In addition, we define a subclass of the above devices, 
we call voltage-controlled memcapacitors, when Eqs. (TTOl 
and ( flTt reduce to 



V c (t) = C 



q(t) = C 



V c (r) dr 



Vc(t) 



q (r) dr 



?(*), 



(15) 



(14) 



where in the above two equations, the lower integration limit 
(initial moment of time) may be selected as —00, or if 
/_ Vc( T )dT = and J_ q(r)dr = 0, respectively. 

Properties - It follows from Eq. (TlOb that the charge is zero 
whenever the voltage is zero. Note, however, that in this case, 
7 = does not imply q — 0, and thus this device can store 
energy. Therefore, the latter can be both added to and removed 
from a memcapacitive system. 

From a microscopic point of view a change in capacitance 
can occur in two ways: i) due to a geometrical change of 
the system (e.g., a variation in its structural shape), or ii) in 
the quantum-mechanical properties of the carriers and bound 
charges of the materials composing the capacitor (manifested, 
e.g., in a history-dependent permittivity e), or both. 

In either case, inelastic (dissipative) effects may be involved 
in changing the capacitance of the system upon application 
of the external control parameter. These dissipative processes 
release energy in the form of heating of the materials com- 
posing the capacitor. However, this heat may not be simply 
accounted for as a (time-dependent) resistance in series with 
the capacitor. 

Similarly, there may be situations in which energy not 
from the control parameter but from sources that control the 
equation of motion for the state variable, Eq. dTTb . is needed to 
vary the capacitance of the system (e.g., in the form of elastic 
energy or via a power source that controls, say, the permittivity 
of the system via a polarization field). This energy can then be 
released in the circuit thus amplifying the current. Therefore, 
Eqs. ( [Tol l and (flTt for memcapacitive systems postulated above 
could, in principle, describe both active and passive devices. 
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However, starting from a fully discharged state, the amount 
of removed energy from a passive memcapacitive system can 
not exceed the amount of previously added energy. Mathemat- 
ically, this property can be written in the form 

t 

U c = J V c {t)I{t)(1t > (16) 

to 

under the condition that at t = to no energy is stored in the 
system. Eq. ( TToT l should be valid for any form of the control 
parameter, such as the voltage Vc(t) applied to the passive 
memcapacitive system. Moreover, in memcapacitive systems 
with dissipative processes (such as heating), the equality sign 
in Eq. (fToT ) for t > to is not realized (assuming Vc{t) ^ 0). 

Fig. |3] shows a schematic memcapacitive system hysteresis 
loop passing through the origin. The shaded areas give the 
energy C/j added to/removed from the system. Note that this 
energy is associated with some internal degree of freedom, 
i.e., related to some physical elastic or inelastic processes that 
accompany the change in conductance. The system is passive 
if U\ + U2 = 0, dissipative if U\ + U2 > and active if 
U x + U 2 < 0. 

Finally, since the state equation (fTTT > has (by assumption) 
only a unique solution at any given time t > to, then if Vc(t) 
is periodic with frequency u>, namely Vc(t) = Vq cos(2-7ra;i), 
then the q — Vc curve is a simple loop passing through the 
origin, namely there may be at most two values of the charge 
q for a given voltage Vc, if we consider a voltage-controlled 
device (see Fig. [4j, or two values of the voltage Vc for a 
given charge q, for a charge-controlled system. This loop is 
also anti-symmetric with respect to the origin if, for the case of 
Eq. mm . C(x,V c ,t) = C(x,-V c ,t) mdf(x,V c ,t) = 
f(x,-V c ,t). 

Like the case of memristive systems, a memcapacitive 
system behaves as a linear capacitor in the limit of infinite 
frequency, and as a non-linear capacitor in the limit of zero 
frequency, assuming Eqs. ( fTTT i and (Tl~3T > admit a steady-state 
solution. The origin of this behavior rests again on the system's 
ability to adjust to a slow change in bias (for low frequencies) 
and the reverse: its inability to respond to extremely high 
frequency oscillations. 

Examples - There are instances in which the capacitance 
C has been found to have a hysteresis-type dependence on 
the applied voltage ll2D . Il22l . ||23l . These cases are related 
to nanoscale capacitors in which interface traps ll2p or em- 
bedded nanocrystals l22l . ||23l are responsible for memory 
effects. Even though characteristic hysteresis properties of 
these systems were not discussed in the framework of the 
memcapacitive systems suggested in this paper, it is clear that 
they belong to the present classification. 

Below, we provide a model example of a memcapaci- 
tive system which could be used to simulate experimental 
realizations of memcapacitors. To be specific, let us con- 
sider a voltage-controlled memcapacitive system according 
to Eqs. ([Tol l and (fTTb . We focus on the electronic circuit 
shown in the inset of Fig. [4] This scheme consists of a 
memcapacitive system C connected to a source of ac-voltage 
V(t) = Vq s'm(2irLut). The capacitance C is chosen to vary in 
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Fig. 4. (Color online) Simulation of a circuit (shown in the inset of (a)) 
composed of a voltage-controlled memcapacitive system and a resistance. 
The following parameter values were used in this calculation: C\ = IfiF, 
C 2 = 20/iF, V T = 0.5V, a = -50mF/(V s), /3 = 2a, Vo = IV. 



the range between its minimum and maximum values (Ci to 
C2) and changes according to the equations 

C = x (17) 
± = - ((3V C + 0.5 (a- 13) [\Vc + V T \ - \V C -V T \]) 
x [9 (Vc) 9(x-d)+e (-Vc) 9 (C 2 - x)] . (18) 

Here, Vc is the voltage drop on the memcapacitive system, 
Vp is a threshold voltage, a is a variation rate for \Vc\ < Vp 
and (3 is a variation rate for \Vc\ > Vr- The interpretation 
of the coefficients a and (3 is the same as in Fig. [2] The 
response of the circuit is simply described by the equation 
V(t) = Vc = q/C where C can be found using Eqs. ( 1171181 1. 
It is important to remember that in modeling of more complex 
circuits the chain rule of differentiation [24] should be used. 

Fig. |4] represents results of our numerical simulations. 
Typical features discussed above such as capacitance hystere- 
sis, low-frequency and high-frequency behavior, and pinched 
hysteresis loop can be clearly distinguished. We also checked 
the energy conservation for the given circuit and found that 
for the present model the memcapacitive system behaves as a 
passive device. 
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IV. Meminductive systems 

Let us now introduce the third class of memory devices. Let 
us define the "flux" 



c/>(t) = / V L (t')dt', 



(19) 



where Vj,(t) is the induced voltage on the inductor (equal to 
minus the electromotive force). 

Definition - We call an nth-order current-controlled memin- 
ductive system one described by the equations 



X 



L(x,I,t)I(t) 
f(x,I,t) 



(20) 
(21) 



where L is the meminductance, and an nth-order flux- 
controlled meminductive system the following 



I(t) = L~ x (x, <f>, t) 4>{t) 



(22) 
(23) 



with L^ 1 the inverse meminductance. 

We call current-controlled meminductors a subclass of the 
above devices, when Eqs. (f20b and (fJTJ reduce to 



4>{t) = L 



I (t) dr 



I(t) 



(24) 



and a flux-controlled meminductor when Eqs. 
can be written as 



and 



I(t) = L- 



<p (t) dr 



(25) 



In the above two equations, the lower integration limit may be 
selected as — oo, or if J_ x 7(r)dr = and /_ ^dr = 
0, respectively. 

Properties - For the sake of definiteness, we shall consider 
the current-controlled meminductive systems. Taking the time 
derivative of both sides of Eq. (|20| > yields 



TA d(j> dl r dL 

VL = Tt= L Tt +I H- 



(26) 



The second term on the right-hand side of the above equation 
is an additional contribution to the induced voltage due to a 
time-dependent L. The energy stored in the current-controlled 
meminductive system can be thus calculated as 



u L {t) 



V L (T)I( T )dT 



dl 

J dt 



r dL 
dt 



I(r)dT. (27) 



in 



When L is constant, we readily obtain the well-known 
expression for the energy Ul = LI 2 /2. Usually, this Ul is 
interpreted as the energy of the magnetic field generated by the 
current. Similarly to the passivity criterion of memcapacitive 
systems (see Eq. ( fTol l) Eq. ( 1271 ) provides the passivity criterium 



of meminductive system if at t — t the latter is in its minimal 
energy state and U^t) > at any time. 

Finally, as in the case of memristive and memcapacitive 
systems, it is expected that when the external control parameter 
is a periodic function of time with frequency lj, namely I(t) = 
Io cos(27rojt), the 4> — I curve is a simple loop passing through 
the origin, namely there may be at most two values of the 
flux <f) for a given current I. This loop is anti-symmetric with 
respect to the origin if, for the case of Eq. (12012 11 1. L (x, I, t) = 
L (x, -I, t) and / (x, I,t) = f (x, -I, t). 

Also, the magnetic permeability can adjust to slow periodic 
variations of the current, while it cannot at high frequencies, 
so that this device behaves as a non-linear inductor at low 
frequencies (assuming Eqs. d2Tb and (f23b admit a steady-state 
solution) and a linear inductor at high frequencies. 

Examples - In electronics, inductors are primarily of a 
solenoid type consisting of a coil of conducting material 
wrapped around a core. The inductance of a solenoid is 
proportional to the relative permeability p r of the material 
within the solenoid and also depends on the geometrical 
parameters of the system. The simplest way to introduce a 
memory effects in such a system is to use the core material 
whose response to the applied magnetic field depends on its 
history. As an example, we can think about ferromagnetic 
materials exhibiting a magnetic hysteresis such as iron. In fact, 
an electronic circuit with such inductor was recently analyzed 
l25l . Another way to introduce memory effects is by varying 
the inductor's shape. 

A circuit model for a meminductive system can be formu- 
lated similarly to the models of memristive and memcapacitive 
systems discussed above. In particular, we can imagine a 
structure whose inductance increases up to a certain value 
when current flows in one direction and decreases down to 
another value if the current flow direction is reversed. This is 
a straightforward calculation which does not add much to the 
analysis of this problem. Therefore, we do not report it here. 

V. Conclusions 

We have extended the notion of memory devices to both 
capacitative and inductive systems. These devices have specific 
properties that appear most strikingly as a pinched hysteretic 
loop in the two constitutive variables that define them: current- 
voltage for memristive systems, charge-voltage for memcapac- 
itive systems, and current-flux for the meminductive systems. 
Many systems belong to these classifications, especially those 
of nanoscale dimensions. Indeed, with advances in the minia- 
turization of devices, these concepts are likely to become more 
relevant since at the nanoscale the dynamical properties of 
electrons and ions may strongly depend on the history of the 
system, at least within certain time scales [16|. 

Apart from the obvious use of these devices in non-volatile 
memories, several applications can be already envisioned for 
these systems, especially in neuromorphic devices to simulate 
learning, adaptive and spontaneous behavior. For instance, the 
identification of memristive behavior in primitive organisms 
such as amoebas |10|, opens up the possibility to relate 
physiological processes that occur in cells with the theory of 
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memory devices presented here. Along similar lines, one could 
envision simple models that identify memory mechanisms in 
neurons and use these memory devices to build such models 
in the laboratory. 

Therefore, due to their versatility (including analog func- 
tionalities) the combined operations of these memory devices 
in electronic circuits is still largely unexplored, and we hope 
our work will motivate experimental and theoretical investiga- 
tions in this direction. 
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